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Abstract 

I review the constraints placed on relativistic pulsar winds by comparing optical 
and X-ray images of the inner Crab Nebula on the one hand with two-dimensional 
MHD simulations on the other. The various proposals in the literature for achieving 
the low magnetisation required at the inner edge of the Nebula, are then discussed, 
emphasising that of dissipation in the striped-wind picture. The possibility of direct 
observation of the wind is examined. Based on the predicted orientation of the po- 
larisation vector, I outline a new argument suggesting that the off-pulse component 
of the optical emission of the Crab pulsar originates in the wind. 
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1 Introduction 



Observations of pulsar nebula e have improved drama tically over the past few 
years . In pa rticular, the X-ray ( Weisskopf et al. . 200fi ) and optical ( Hester et al 



199i l2002h images of the Crab Nebula have revealed intriguing insights into 
the dynamics of this object that have motivated detailed theoretical studies. 
These, in turn, have sharpened the constraints on the possible physical con- 
ditions at the inner boundary of the Nebula, which encloses both the pulsar 
and an apparently empty region surrounding it. This paper concentrates on 
the physics of the apparently empty region, that I will call the pulsar wind, 
because it contains out-flowing material, moving, presumably, at supersonic 
velocity. How ever, it also contains large amplitude waves of the s ame period 



as the pulsar ([Rees and Gunnl . Il974t iKundt anc 



stantial dc component of magnetic field ([Pacini 



Kxotscheck , 19801 ) and a sub- 
1968lh 



Whereas the outer boundary of a pulsar wind is, at least in some cases, di- 
rectly accessible to observation, constraints on the physical conditions where 
the wind is launched are difficult to find. A common scenario is that the pulsar 



Preprint submitted to Elsevier Science 



2 February 2008 



magnetosphere drives an almost radial wind that becomes supersonic at some 
point not too far outside the light cylinder and that contains a current sheet 
separating the two magnetic hemispheres defined by the dipole component of 
the field at the stellar surface. However, althoug h this is an att ractive model 
to work on, alternative pictures are possible because our 



knowledge of the magnetosphere is incomplete. Even the term magnetosphere 
itself can cause confusion. It is generally used to identify a region in which 
external electromagnetic fields dominate the dynamics of a plasma. For pul- 
sars, this certainly applies inside the light cylinder. But the definition becomes 
ambiguous further out, where the magnetic field weakens and the outflow ac- 
celerates. This is because electromagnetic fields may still dominate the energy 
flux even in a supermagnetosonic flow, in which the inertia of the plasma is 
important. Here, I use the term magnetosphere to refer to regions close to the 
neutron star where the bulk radial flow is still submagnetosonic, including, of 
course, those parts within the pulsar light-cylinder that corotate with the star. 
Understood in this sense, it is generally thought that the wind is dark, but 
the magnetosphere is responsible for the pulsed emission detected, in the case 
of the Crab, from the radio to gamma-ray bands as a point source apparently 
located at the position of the neutron star. As discussed in Section 5 this is 
not necessarily the case. 

In this paper I will briefly review the constraints that can be placed on the 
physics of the pulsar wind by MHD simulations of the Crab Nebula. I will 
then consider current theoretical ideas on how these constraints could be met, 
and finally look at the prospects for observing the wind directly. 



2 MHD simulations 



Axisymmetric, relativistic MHD simulations of the wind and the b ubble it in- 



flates downs trea m of the termination sho ck have been performed bv lKomissarov and Lvubarskv 
f)200.il2004h andEeT Zanna et alJ (J200J). The results are broadly similar. Pro- 
vided the energy flux in the wind is assumed to be concentrated towards 
the equatorial plane, the termination shock is highl y oblate - - an a spect 
that had been deduced using analyt ical arguments by Lvubarskv ( 2002j ) and 
Bogovalov and Khangoulian (|2002al lbl). If, in addition, the magnetisation pa- 



rameter a = B 2 / (4irw) (here B is the magnetic field in the plasma rest frame 
and w the proper enthalpy density) is tuned to the correct value, th en the mag- 
netic pinch in the shocked plasma collimates the flow, as predicted flLvubarskv , 
2002i ) . and a torus-plus-jet like structure similar in appearance to that seen in 



the optical and X-ray images of the Crab Nebula is reproduced. The anisotropy 
of the energy flux and the tuning of the magnetisation are critical; the simu- 
lations are less sensitive to other parameters such as the Lorentz factor of the 
flow and the angular distribution of the particle flux. 



2 



The anisotropy of the energy flux should follow from the boundary conditions 
on the magnetic field at the stellar surface. But our knowledge of this con- 
nection is s canty. The on ly known exact solution is that of the force-free split 
monopole (|Michell . ll97i . In this case, spherical symmetry of the magnetic 
field at small radius yields an energy flow per solid angle, dL/dfl, that is con- 
centrated into the equatorial plane, with a dependence on colatitude 9 given 
by 



The poloidal field in this solution is purely radial, and there are no regions 
of closed field. Nevertheless, a non-axisymmetric version c an be cons t ructe d 
by simply inclining the plane used to "split" the monopole (jBogovalov . 1999). 
In a realistic case, it is expected that the closed field line regi ons fill much of the 
magne t osphere ins i de the light cylinder [see the discus sion in lContopoulos et al 
f)l999h . lUzdenskvl (|2003h and lOoodwin et all (|2004h ]. Consequently, there is 
no reason to suppose that the angular distribution given by Eq. (1) is a good 
approximation. Nevertheless, all simulations published so far start off with 
essentially this distribution. 



There is less consensus concerning the magnetisation parameter a. This quan- 
tity is, in principle, also a function of angle, since it describes the way in which 
the energy flux at a particular colatitude is divided between Poynting flux and 
particle-born energy. In the exact force-free solution a is formally infinite. But 
a simulation with large a would fail to look anything like the Crab Nebula. In 
fact, for a greater than about 1% in the equatorial region, the jet appears too 
fast and strong and the equatorial outflow is suppressed. For smaller a, some of 
the detail ed properties of the outflow match obse rved features such as the cen- 
tral k not ( Hester et al. . 19951 ) quite convincingly ( Komissarov and Lvubarskv . 
2004). Thus, it appears that the magnetisation parameter must be very small 
in the equatorial wind. This is a well-known and puzzling property, as discussed 
below. It is sometimes termed the "cr paradox", and w as already recognised 
in sp herically symmetric models of the Crab Nebula ((Kennel and Coroniti . 
1984). On the other hand, if a is everywhere small, the two-dimensional simu- 
l ations resemble the purely hydrody namical case, resulting in a "lonely torus" 
(IKomissarov arid Lvuharskvl. 120041). 



The value of a in the energetically less important, polar parts of the outflow is 



uncertain. 



Del Zanna et al.l ((20041 ) concentrate on the case in which a is inde- 



pendent of angle, with a current sheet in the equa torial plane, as in the case of 
the "aligned " split monopole. On the other hand, IKomissarov and Lvubarskv 
f)2003Ll2004h are motivated by the scenario where low magnetisation is achieved 
by complete dissipation of the oscillating component of the magnetic field in 
the "striped wind" of the oblique split-monopole solution ((Coroniti . 1990; 
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Fig. 1. The striped pulsar wind. A magnetic dipole embedded in the star at an 
oblique angle to the rotation axis introduces field lines of both polarities into the 
equatorial plane. The figure shows the intersection of the current sheet separating 
these regions and the equatorial plane. In the inset, an almost planar portion of this 
sheet (dashed line) is shown, together with the magnetic field lines, assuming they 
undergo reconnection. 



z 




Table 1 

Dependence of the bulk Lorentz factor T of the pulsar wind on radius for different 
dissipation mechanisms. r max is the radius at which dissipation is complete, j"l is the 
radius of the light-cylinder, fi = L/(N±mc 2 ) is Michel's parameter (with L the spin- 
down luminosity, and N± the pair loss rate in the wind) and L = Lir 3 e 2 /(m 2 c 5 Q, w ) 
is a dimensionless measure of the spin-down luminosity divided by the solid angle 
O w occupied by the wind (L 1.5 x 10 22 for the Crab pulsar). L is approximately 
equal to the square of the pole to equator potential in units of the electron rest-mass. 

Slow dissipation Tearing-mode Fast dissipation 

T oc r 1 / 2 r cx r 5 / 12 T oc r 1 / 3 

r max /r L = L 1 / 2 r mSuX /r L = ^ 4 / 5 L 3 / 10 r max /r L = ^ 2 



Lvubarskv and Kirkl . 120011 ). On the equator, the dc component of the field in 
this solution vanishes for all obliquities, and so o = here after dissipation of 
the stripes. But, for all except the exactly perpendicular case, the stripes are 
confined to a band of latitudes around the equator, implying a high value of 
cr close to the poles. Interestingly, both groups agree that the average value 
of a, i.e., the ratio of the total luminosity in Poynting flux to the total in 
particle-born flux, should be a few percent. 
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3 The a paradox 



The exact solution with a — > oo found by iMichell (|l973h has radial poloidal 
field. Its large a counterpart also appears neithe r to collimate significantly, 
nor to accelerate ( Bogovalov and Tsinganosl . 1999). Thus, if, as expected, a is 
large close to the star, it should continue to be so out to the termination shock. 
The small magnetisation implied by observations completes the a paradox. 

There are several possible escape routes: 

(1) Before the adve nt of the high-resolution images of the Crab Nebula, 
iBegelmanl (|l998l ) suggested that a remains large even outside the ter- 
mination shock, and that the magnetic energy dissipates in the Nebula. 
However, this now seems to be ruled out by the good agreement of low a, 
ideal MHD simulat ions and observations in the optical and X-ray bands. 

(2) iLvubarskvi (|2003bh also proposed that a remains large in the wind, but 
that the energy carried in the oscillating component of the magnetic 
field dissipates in a transition region that can be considered as a special 
kind of termination shock. Such a configuration leads to the same initial 
conditions for the MHD simulations as does a low a wind, provided the 
transition region remains thin. This is an interesting suggestion that needs 
further investigation. In particular, it would be important to identify any 
spectral signature that could depend on the position of the termination 
shock. 

Some collimation and, therefore, acceleration and conversion of magnetic 
into particle-born energy might be possi ble if the initi a l dist r ibution of 
poloi dal flux is sufficiently anisotropic (|Chiueh et al. . 19981 : Iviahakisl . 



(3) 



20041 ). However, whether this is a realistic possibility can only be de- 



termined by co nstru cting a model of the inner magnetosphere. 
(4) ICoroniti ( 1990) and Michell (1994) proposed that the oscillating compo- 
nent of the magnetic field in the striped wind picture should gradually 
damp by dissipation of the magnetic energy at the field line reversals. 
However, these early papers did not include the associated acceleration 
of the wind and the resulting dilation of the dissipation timescale. 

Routes (2) and (3) await further theoretical development, and I now turn to 
a discussion of route (4) . 



4 Dissipation mechanisms 



To assess the importance of wave damping as a possible solution of the a prob- 
lem it is first necessary to identify the nature of the wave. Close to the pulsar, 
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electr omagnetic modes do not propag ate (|Usovl . 1197,4 iMelatos and Melrosd . 
1996). On the other hand, both the fast magnetoson ic wave ([Lvubarskvl . 
2003cJ ) and the entropy wave ([Lvubarskv and Kirkl . 120011 1 are possible. In prin- 
ciple, t hese may conv ert spontaneously into electromagnetic modes at larger 
radius ( Melatosl . 1998| ). but the conditions under which this happens have not 
yet been fully inves t igated . If the waves persist as subluminal modes, damping 
is inevitable ( Usov . 19751 ). 



The rate of damping of the striped wind pattern and, therefore, the viability 
of escape route (4), depends on the speed with which magnetic reconnection 
can proceed in a relativistic current sheet. These structures (see Fig. 1) dif- 
fer from their more familiar counterparts in solar and magnet ospheric physics 
because they do not permit evacuation of the heated plasm a from the sheet 
as in a Sweet-Parker or Petschek-type reconnection model ((Kirk . 2004). In- 
stead, either the sheet steadily thickens, or plasma is ejected along the induced 
electric field (the ^-direction in Fig. 1). As dissipation proceeds, the current 
sheet separating regions of opposite magnetic polarity fills with hot plasma. If 
this process proceeds slowly (not faster than the timescale needed for a fluid 
element in the wind to double its radius) a small- wavelength approximation 
can be used to analyse the evolution of the stripes. On the other hand, rapid 
dissipation cannot be ruled out, and is difficult to quantify. 



Concentrating on the small-wavelength approximation, the physical picture is 
one in which the sheet remains in approximate pressure equilibrium with the 
surrounding magnetic field, which is dragged into it and annihilated. Spherical 
symmetry must be assumed to make the analysis tractable, but since the flow 
is highly relativistic, pressure imbalance in the 9 direction is unlikely to be 
important, and we can assume the solution holds independently along each 
radius vector. A lower limit on the rate at which annihilation can proceed is 
found by arguing that the sheet width c annot be smaller than the g yro radius 
of a "hot" particle in the confining field ((Lvubarskv and Kirku2001[ ). However, 
this rate is too slow to be of interest in the case of the Crab. A more realistic 
estimate of the annihilation rate follows by assuming it eq uals the timescale 
on which the relativistic te aring-mode instability operates ( Lvubarskii . 1996; 
Kirk and Skiaeraase But, again, this is not fast enough to resolve the 

er paradox for the Crab. An upper limit, subject to the assumptions inherent 
in this method, is given by demanding that the hot sheet expand at most 
sonically. This does provide a viable escape route for the case of the Crab, but 
only if the number of electron/positron pairs carried from the pulsar by the 
wind is significantly larger (> 3 x lO^s -1 ) than the rate fo und in standard 
pair-production calculations (|Hibschman and Aronsl . l2001al fbl). Independent 
support for injection of pairs into the Nebula at this high r ate is provided 



by m odels of the radio to hard X-ray synchrotron spectrum ([Gallant et al 



2002t ). so that we may still have much to learn about where and how a pulsar 
produces electron/positron pairs. 
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For each of these proposed dissipation processes, the acceleration of the wind 
can be described by a similarity solution in which the bulk Lorentz factor 
is proportional to a power of the radius. These results are summarised in 
Table 1. In principle, they make possible a prediction of both the average a 
and the Lorentz factor of a pulsar wind as a function of radius. Given the 
inclination angle of the magnetic dipole axis, predictions can also be made of 
the dependence of these quantities on colatitude. 



5 Observing the wind 



In the case of gamma-rays, an observational upper limit on the wind emis- 
sion yields a firm constraint on the pulsar wind. Even if the plasma remains 
completely cold, if its bulk Lorentz factor is large, it will radiate by Comp- 
ton upscattering the ambient target photons. Using thermal X-rays fro m the 
surface of the Crab pulsar as targets, Bogovalov and Aharonian ( 20001 ) inter- 



preted the detected unpulsed flux in TeV gamma-rays as an upper limit on 
the contribution of the wind, and derived a lower limit of 5?"l on the radius 
at which the Lorentz factor of the Crab wind can reach 10 6 . This constraint 
is satisfied by all the models presented in Table 1. 

In the optical and X-ray images the region around the Crab pulsar appears 
dark out to a distance of about 12 arcsec. However, this does not mean the 
wind does not radiate, but merely that along a line of sight that passes between 
about 0.5 arcsec and 12 arcsec of the pulsar, the emission of the wind is beamed 
away from the observer. This occurs quite naturally in a relativistic, radial 
wind. If the bulk Lorentz factor is proportional to r q with q < 1 — as in the 
cases considered in Table 1 — the most stringent constraint is imposed at the 
termination shock. For the emission to be beamed away from the observer at 
that radius requires T > 12 arcsec/0.5 arcsec = 24. This is much smaller than 
the Lorentz factor usually assumed for the plasma entering the termination 
shock. Consequently, even if it radiates at all radii, the wind should appear to 
contribute only to the point source, leaving the surrounding region dark. 

An analogous argument can be advanced concerning pulsation of the wind 
emission. If the striped pattern gives rise to a phase dependence of the emis- 
sivity, for example, to increased emission at the field reversals, the radiation 
from radiu s r will appear pu lsed at the rotation period of the star provided 



from radiu s r will appear pu 
r/r L < T 2 (|Kirk et all [2002). 



Has this radiation been observed? To answer this question, one must be able to 
untangle the contribution of the stellar surface, magnetosphere and wind, all of 
which are co-spatial (and point-like) when observed with current instruments. 
This is a difficult task, especially given the uncertainties associated with mod- 
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elling the pulsed emission. However, one characteristic of synchrotron emission 
from the wind that does not appear to be shared by any magnetospheric model 



2P 

fjDvks et all 12004 iKaspi et all I2004T ) is the direction of linear polarisation of 
the off-pulse or "dc" component. In the striped wind, the polarisation vec- 
tor of radiation emitted close to the field reversals is not well constrained, 
but outside of these potentially pulse-producing regions, the magnetic field 
is purely toroidal. This predicts that the electric vector of the dc component 
of synchrotron radiation should lie along the projection onto the sky of the 
rotation axis of the neutron star. In the case of the Crab pulsa r, the optical 



emission appears to have precisely this property ( Kellner . 2002j ). 
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